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The magnetic excitations of the double exchange (DE) model are usually discussed in terms of an 
equivalent ferromagnetic Heisenberg model. We argue that this equivalence is valid only at a quasi- 
classical level — both quantum and thermal corrections to the magnetic properties of DE model differ 
from any effective Heisenberg model because its spin excitations interact only indirectly, through the 
exchange of charge fluctuations. To demonstrate this, we perform a novel large S expansion for the 
coupled spin and charge degrees of freedom of the DE model, aimed at projecting out all electrons 
not locally aligned with core spins. We generalized the Holstein-Primakoff transformation to the 
case when the length of the spin is by itself an operator, and explicitly constructed new fermionic 
and bosonic operators to fourth order in 1/y/S. This procedure removes all spin variables from the 
Hund coupling term, and yields an effective Hamiltonian with an overall scale of electron hopping, 
for which we evaluate corrections to the magnetic and electronic properties in 1/S expansion to 
order 0(1/5'^). We also consider the effect of a direct superexchange antiferromagnetic interaction 
between core spins. We find that the competition between ferromagnetic double exchange and an 
antiferromagnetic superexchange provides a new example of an "order from disorder" phenomenon 
— when the two interactions are of comparable strength, an intermediate spin configuration (either 
a canted or a spiral state) is selected by quantum and/or thermal fluctuations. 



I. INTRODUCTION 

Many msufimetic systems ofpgreat experimental interest, for example the colossal magnetoresistance (CMR) 
Manganiteall'B and PyroclhoresB, comprise a band of itinerant electrons interacting with an ordered array of localized 
magnetic moments with spin S. In many cases these systems can be modeled as a single tight band of electrons 
interacting with localized core spins by a Hund rule exchange interaction, and can be described by a Hamiltonian of 
the form : 

n = nn + ni+n2 (i) 

where Hq is the bare Hamiltonian for itinerant electrons 

-Ho^-tY, cLcja, (2) 

{ij)a 

Til describes the first Hund rule (Kondo type) coupling between localized spins and itinerant electrons 

"Hi = -Jh ^ St.Si = ^ c\^(^2(iCti3 (3) 

i Q/3 

and TL2 describes a superexchange between localized spins 

H2=-J2^5,5, (4) 

Here the indices i and j run over lattice sites, a and [3 over electron spin states, and components of cjap are Pauli 
matrices. 

Our primary interest will be to study the limit in which the Hund's rule coupling Jh is positive and is much larger 
than both J2 and t. In this case itinerant conduction electrons must be locally aligned with the core spin on any site. 
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As the kinetic energy of electrons is minimal when all electrons are parallel, the core spins are also all parallel, i. e. the 
ground state is a ferromagnet. The ferromagnetic interaction between core spins mediated by conduction electrons is 
often referred to as "double exchange", and the model described by Ho + 'Hi, in the limit t/Jn 0, is often referred 
to as the double exchange ferromagnet (DEFM). 

Many authors have demonstrated that for classical spins, the fermion-mediated ferromagnetism is described by the 
effective nearest-neighbor Heisenberg model 

-Ji^f,.f, (5) 

(iJ) 

with T = S + x/2, where x is the electron density, and ferromagnetic exchange integral 

J. - ^ (6) 

where t is the kinetic energy per bond in the lattice, which is a function of the electron doping x. In the limit of low 
electron density x — > 0, i ~ tx. 

The aims of this paper are two-fold. Firstly, we discuss the extent to which the DEFM is equivalent to the nearest- 
neighbor Heisenberg ferromagnet. We argue that this equivalence holds only at infinite 5*, when spin waves are 
non-interacting quasiparticles. We demonstrate that the interaction between spin waves in the DEFM is qualitatively 
different from that in the Heisenberg ferromagnet, for which spin waves are exact eigenstates of the Hamiltonian. 
This difference gives rise to different forms of quantum and thermal corrections to the spin wave dispersion. In 
particular, we show that spin excitations in the DEFM have a finite lifetime even at T = 0, i. e. that they cannot be 
true eigenstates of the Hamiltonian. Physically, the difference between the two models stems from the fact that in a 
DEFM, the dynamics of the bosonic spin wave modes are completely determined by those of the itinerant electrons. 
Just as the existence of a finite density of charge carriers generates a dispersion for spin waves about the groundstate 
of the DEFM, so fluctuations of charge density generate a retarded interaction between these spin waves. Since this 
interaction is proportional to the charge susceptibility of the itinerant electrons, both its value and the dependence 
on momentum and frequency are very different from that in the Heisenberg model. 

Our second goal is to analyze what happens when the ferromagnetic DE interaction competes with the antiferromag- 
netic superexchange. The trade-off between fermion-mediated and direct exchanges caiL in principle, lead to many 
different ground states, and indeed numerical studies suggest a very rich phase diagrama. Here we address the issue 
how the system evolvesJrom a DEFM to and AFM with increasing J2, provided that there is no phase separation. As 
observed by de GennesQ, for classical spins, the configurations which interpolate between FM and AF order are the 
ones in which the neighbouring lattice spins arc misaligned by an angle 9 such that cos 6*72 = J1/J2. This criterium 
can be satisfied by canting the spins into a two-sublattice state shown in Fig. [|a or into a spiral state shown in Fig. 
However, these two configuration are not the only possible as at the classical level we may take any spin of, say, B 
sublattice of the canted phase and rotate it about the direction of magnetization of the A sublattice without altering 
the angle between it and the neighbouring spins. This rotation introduces an infinite set of classically degenerate 
intermediate configurations which interpolate between canted and spiral states. We performed a spin wave analysis 
of the canted phase using the transformation to bosons described below, and indeed found that a local degeneracy 
yields a branch of spin wave excitations with ujsf{q) = for all q. Since it costs no energy to make-an excitation, the 
system cannot distinguish between different states, and is magnetically disordered even at T = OQ'EI. 

This argument, however, does not hold for quantum spins, and we expect that quantum fluctuations enable the 
system to choose its true groundstate. A widely studied example^of this "order from disorder" behavior is provided 
by the highly frustrated 2D antiferromagnet on a Kagome latticeQ. The question remains as to what kind of order is 
preferred in our case. 

To address this issue we analyzed what is the momentum of the instability of a FM configuration. For classical 
spins, as discussed above, DE ferromagnetism can be described by the effective exchange model Eq. (^). In this 
situation, the excitation spectrum in a FM phase is simply 

ujsw{q) = zS{Ji - J2)(l - 7g) (7) 

where z is the lattice coordination number, and 



z 
s 



r ^ — ^ 



where the sum on {6} runs over nearest neighbors vectors. This excitation spectrum vanishes identically for Ji = J2, 
in agreement with the infinite degeneracy at Ji < J2. Suppose now that this degeneracy is lifted by quantum 
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FIG. 1. Two sublattice canted (a) and spiral (b) spin configurations. 



fluctuations, and the spin wave spectrum first becomes unstable at some momentum Q. One can easily demonstrate 
(see below) that if Q = 0, the resulting state is a spiral, ii Q = (tt, tt, tt), the instability leads to a canted state, and if 
Q is in between the two limits, the system chooses some intermediate spin configuration. 

We show that for both small and large electronic densities (x ^ 1 or 1 — .x <C 1), the DEFM becomes unstable 
(with increasing J2) against the two-sublattice canted structure. At intermediate densities (0.05 < a; < 0.68 in 3D), 
the first instability is against a spiral spin configuration. Similar results hold for the 2D case. We also show that 
thermal corrections to ujsw{q) are very different from those in the Heisenberg ferromagnet, and in 3D can be even of 
a different sign at the lowest densities of carriers. For realistic densities, we found that the functional form of the 
thermal correction is approximately the same as in the Heisenberg model, but the overall amplitude is substantially 
increased. 

To proceed with the spin wave calculations, we need a "bosonization" scheme which treats both core spins and 
itinerant electrons on equal footing. Our approach is to transcribe the model Hamiltonian of Eq. (|l]) in terms of 
the "natural" collective coordinates of the Kondo coupling term under the assumption that the size of the core spin 
S ^ 1. This procedure has the advantage of making a clean distinction between different types of excitations (spin 
and charge modes) and of showing clearly the separation of different physical energy scales in the atomic limit. The 
transformation and its derivation are presented in the next Sec. ^ In Sec. |l|, we apply this transformation to the 
DE model and discuss the form of the thermal and quantum corrections to the spin wave spectrum. In this section, 
we also present the results for the damping rate of spin excitations. In Sec. we discuss how quantum and/or 
thermal fluctuations select an intermediate configuration which interpolates between ferromagnetic state at Ji > J2 
and antifcrromagnetic state at vanishing J1/J2. We also discuss here a peculiar re-entrant transition between spiral 
and canted states, which is due to a competition between the selection by thermal and quantum fiuctuations. Finally, 
in Sec. Mwe present our conclusions. 



Since we wish to work in the limit where the Hund's rule coupling Jh ^ t is the largest energy scale in the problem, 
it makes sense to treat the onsite Kondo-coupling exactly, and to introduce hopping between sites a perturbation. 
In the limit where Jh 00, we can go one step further and project out all electrons not locally aligned with the 
core spins. We should also define spin wave excitations such that they are the true Goldstone modes of the order 
parameter, i. e. transverse fluctuations of the composite spin 



We can accomplish these goals by first introducing new Fermi operators fj and pj which create local states with high 
{T = S + 1/2) and low {T = S — 1/2) total spin T, and then generalizing the Holstein-Primakoff transformation to 



II. DERIVATION OF TRANSFORMATION 




(9) 



Q/3 
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the case where the length of the spin T is itself aruoperator, introducing a corresponding bosonic operator a^. Here 
we outline this procedure, which was introduced ina. A comparison with an exact solution on two sites are given inLI. 

Let us start by considering the "Kondo Atom" of a single localized spin and itinerant (spin degenerate) electron 
orbital 



T-Li = - JS s 

We can diagonalizc this Hamiltonian in by introducing the total spin operator 



(10) 



(11) 



such that 



Hi 



J 
2 
J 



= --[T{T+l)-S{S+l)-s{s + l)] 



(12) 



Here the quantum number for the "itinerant" electron s can take on values s — {0, 1/2}, and the total spin quantum 
number T = {S -\- 1/2, S,S — 1/2}, so we find three corresponding energy eigenvalues 



, , JS ^ JS 1 



(13) 



where the state with total spin S and eigenvalue zero is doubly degenerate and corresponds to zero or double occupancy 
of the itinerant electron state. 

We can now define a new quantization axis for a composite spin T and introduce the corresponding "up" (/) and 
"down" {p) fermionic states. The magnitude of T is then 



T = S 



Pf - 



(14) 



This procedure corresponds to making an SU (2) rotation of the electron spin coordinates such that the new quanti- 
zation axis follows the direction of the instantaneous local magnetization. Accordingly, one approach to deriving an 
effective action for the DE model is to work in the SU (2) rotated local frame and introduce a gauge field to describe 
the rotation of coordinates for an electron hopping along any given bond of the lattice (see e. g. Efl). 

For our purposes, it is however more advantageous to use a somewhat different procedure. First we note that to 
leading order in l/S*, the states with T = S'± 1/2 correspond to aligning or antialigning the spin of a single itinerant 
electron with the localized spin on that site. At this level, we can replace the Kondo coupling with a local magnetic 
field and reduce Tii to Tii = Tii + const where 



(15) 



Second, we observe that we can obtain the correct eigenvalues of the full quantum mechanical problem by keeping 
( |r^ ) and modifying the simple "Zeeman-splitting" form of Eq. (|l^) to 



'^^ 



s 



One can easily make sure that Eq. (16) yields correct eigenstates (|l^ ) for any S. The new local interaction 



- {Pfp'p 



(16) 



(17) 



comes from the fact that the eigenvalues of Hi for T = S" ± 1/2 are not symmetric about zero. Physically, it implies 
that if the system has both / and p electrons it can save energy for positive J by putting them on the same site to 
male a total spin of length T = S. 

We see that if J is the largest parameter in the problem and S is large (as we assume) , the fermionic states created 
by the p— operator are separated in energy by J{S + 1/2) from fermionic states created by /—operator, and can be 
safely eliminated from the low-energy sector. This is the obvious advantage of introducing /— and p— operators. The 
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problem which we now have to solve is how to relate these two operators, which describe fermionic states of a composite 
spin T to the "original" up/down fermionic operators C| = u, C| = d introduced with respect to a quantization axis 
of the core spin S. We solve this last problem order by order in l/S* by introducing two sets of Holstein-Primakoff 
bosons — one a for a core spin S and another a for the composite spin T. Introducing the quantization axis z for the 
core spin and the corresponding fermionic states C| = u, C| = d, we have 



5^ 
S+ 



S — a* a 



25\ 1 



ai a 
25 



a 



S- = V2Sa\ I 



(18) 
(19) 

(20) 



and 



T+ = S- 



2 

+ u^d 



= T — a^a 



2T\ 1 



2T 



2T 



(21) 
(22) 

(23) 



where [a, 



[a, d^^l 



1. The relation between a, u,d and d, / and p then follows from Eq. (n4h on the magnitude 



of the composite spin, commutation relation for fermionic and bosonic operators and the constraints that the total 
number of particles be conserved 



v)u - 



dU 



(24) 



and that double or zero occupancy of the electron orbital correspond to the same two sets of states state in either 
representation 



fp 



etc. 



(25) 



This set of constraints is sufficient to uniquely define the transformation, and after lengthly algebra we find : 



f = u 
p — d 



1 

1 - 
1 - 



dU- 



45* 
a + d 

Is 



45* 



Zu'^u-bdU+2u'^udW\ u^d 
32^2 J 



32S'2 

— 2uu^a^ a — a^aa 



1 

" Is 

da'' 



d^i 



2S 45 



32S'2 



2S 

ua 



'-AS 



'-4S 



(26) 
(27) 
(28) 



To verify the transformation, we checked that it correctly reproduces the Kondo coupling term on a single site. In 
terms of original a, u, and d operators, the Hamiltonian for a Kondo "atom" reads : 



Hi 



JS 



—ad^u 




(29) 



Substituting the transformation, we indeed found that in terms of new d, / and p it exactly reproduces Eq. (|16[). 

For purposes of developing a calculation scheme based on these coordinates we need the corresponding inverse 
transformation. This is found to be identical up to the sign of terms at odd order in 1/VS : 



d = p 



1 - 



45 
d^d + p^p 

45 
d^d + //t 

45 



f^P 



3252 

3p^p + 3d^d — 2p^pd^d — d^d^dd 
3252 

3//t + 3dtd - 2//tdtd - d^dtdd 



^-45 
pd^ 



P^f dd 



2545 



3252 



25 
/5 



^-45 
1 

~ 45 



(30) 
(31) 
(32) 
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The inverse transformation Eqs. (30-32|) can be applied whenever the Hamiltonian Eq. has a magnetically ordered 
ground state, and since all Fermi and Bose operators are well defined, is an ideal starting point for constructing 
a diagrammatic perturbation theory of spin and charge excitations in such a system, and for calculating response 
functions of the DE model in a controlled way. l-l 
For the simplest possible case of two core spins sharing a single electron, the model Eq. ([l]) can be solved exactlyn, 
and one of usjaas checked explicitly that our expansion scheme reproduces all features of the exact solution, order by 
order in 1 / ^/Su. 



III. SPIN WAVE EXCITATIONS IN THE DOUBLE EXCHANGE MODEL 



We now return to the DE model and analyze the form of the spin wave dispersion using the transformation derived 
above. As discussed, for a ferromagnetic Kondo coupling Jh and t/Jn — > 0, p— operators can be safely neglected. 
The effective Hamiltonian for the DEFM is therefore found by substituting the inverse transformation Eqs. (j3C|-[3^) 
into Eq. (|l|) and dropping all p operators. Corrections at finite t/Jn will be discussed elsewhere. Substituting the 
transformation into the DE Hamiltonian Ti. = Ti.Q + Tii and expanding in 1/5 we obtain in terms of the new / and d 
operators 



n = 

Tin — 



n'Q + H2 + ni + o{i/s^) 



7^4 



zt 
32S'2 



On Fourier transform we obtain 



fci 



^2 = E "24/lV2a3«4'5l+3-2-4 



1 



7^4 = 



4(5' 
1 1 



1 + ^ ) (ei+3 + £2+4) - (ei + €2) 



3252 iV2 



3+4-2-5-6 



ki...ke 



^256 — ^ [fl+3-5 + £1+3-6 + £1+4-5 + £1+4-6 



4ei 



£2+6-4 + £2+6-3 + £2+5-4 + £2+5-3 ^ 462] 



(33) 
(34) 

(35) 
(36) 

(37) 

(38) 

(39) 
(40) 

(41) 



All energy scales are set by the electron dispersion, which for the simple tight binding kinetic energy term Eq. ^ , is 
given by ek = -ztjk- 

We see that the original problem of spinfuU electrons interacting with quantum mechanical localized spins reduces 
to a single band of spinless fermions interacting with a reservoir of (initially dispersionless) bosonic spin modes. For 
this last problem we can evaluate all quantities of physical interest diagrammatically, starting from the bare bosonic 
Green's function given by 



and the bare / electron Green's function 



C/o(fc, iiOn) = {iuJn - £fc + A*) ^ 



(42) 



(43) 



We note that since operators and d describe fluctuations of the composite spin, the pole in the fully renormalized 
bosonic propagator r>(g, ifln) coincides with the pole in the transverse spin susceptibility, i. e. X>(g, ifln) describes true 
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a) 



c) 




d) 




FIG. 2. a)-d) Diagrams contributing to spin wave self energy in Ferromagnet to order 0{1/S^). Only diagrams a) and b) 
are physically relevant. Diagram c) cancels unphysical leftovers from the diagrams a) and b), and diagram d) accounts for 
thermal renormalization of the chemical potential and the exchange integral Ji but does not change the form of the spin wave 
dispersion. 



spin wave excitations. Higher-order terms in the Holstcin-Primakoff expansion of the composite spin T in terms of 
and S, only give rise to the incoherent background in the spin susceptibility but do not affect the pole. This separation 
between the pole and the incoherent background only makes sense if the damping of spin waves is negligible small. 
We will show that the spin wave damping appears only at 0(1/5''^) such that to order 0(1/5^) (to which we will 
perform controlled calculations) , nonlinear terms in the Holstein-Primakoff transformation for T can be neglectedt2l. 

The dispersionless form oi'D{q^i^n) does not survive self-energy corrections arising from interaction with Fermions 
— the bosonic self-energy depends on momentum g, and this dependence gives rise to a dispersion of the spin wave 
excitations. Physically, this dispersion is generated by the fact that any departure from perfect ferromagnetic ordering 
of composite spins costs kinetic energy. The form of this dispersion should be appropriate to a ferromagnet on a lattice, 
i. e. it should have a set of Goldstone modes with energy uJsw{q) scaling as in the zone center, and be continuous 
across the zone boundary. 

The perturbation theory for the bosonic propagator D{q, 51) is straightforward. We have D~^{q, ft) — il— ^(q, O). 
The lowest-order (in 1/S) contribution to S(q, f2) comes from a single loop of fermions, and this evaluates to 



where 



E(i'^)(g) = zJi5[l-7(«)] 



(44) 



(45) 



The result for Ji can be rewritten as Ji — t/ (25*^) where t is the magnitude of the Fermionic kinetic energy per bond 
on the lattice. 

The classical spin wave dispersion of the DEFM, Eq. 44 corresponds exactly to what we would expect for a nearest 
neighbour Heisenberg model with exchange integral Ji . The tsauivalence of the DEFM and the nearest neighbour 
Heisenberg model at 0(1/5) was first discussedJay de Genneg3 and results for classical spin wave dispersion have 
been rederived by many other authors sincetiHi3ii3. 

We can translate the energy scale for zero temperature spin wave dispersion into a mean field transition temperature 
for a Heisenberg FM using the relation 



zS{S+l) 

-tc — Z "^1 



(46) 



giving approximate transition temperature for a d dimensional quarter filled cubic lattice with bandwidth These 
estimates are in surprisingly good correspondence with transition temperatures for the real CMR materials. 
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FIG. 3. The representation of the diagram Fig gb as an effective four-boson interaction mediated by the charge suscpetibihty 
(particle-hole polarization bubble). 



We now demonstrate that the DE and Heisenberg models are not equivalent beyond 0{1/S). To check this, we need 
to go a step further and compute the spin wave dispersion to order 0(1/5^). For this, we must include the 0(1/5*^) 
contribution of the one-loop diagram Fig. ^ and also two-loop self-energy diagrams. The two-loop diagrams are 
presented in Fig. |^-d. The first of these diagrams represents the bosonic self-energy due to the effective four-boson 
interaction mediated by the Pauli susceptibility of ferniions. The diagram in Fig. 0c is the first-order contribution 
from the six-fold, l/S*^ term in Eq. (p7[). This diagram is physically uninteresting, and is only necessary to restore 
the Goldstone theorem. The diagram Fig. ^d accounts for thermal renormalization of the chemical potential and the 
exchange integral Ji but does not change the functional form of the spin wave dispersion. We neglect this diagram 
below. 

We see that from physics perspective, the relevant diagram is Fig. ^d. This diagram can be thought of as a first 
order bosonic self-energy due to effective four-boson interaction mediated by fermions (see Fig. ||a-b). The relation 
between DEFM and the Heisenberg model then can be readily understood-an general grounds. Indeed, in a HEM, the 
four-boson vertex does not depend on frequency and scales as Jiq^P (Ref.EZI) where q,l 1 are the bosonic momenta. 
In the DEFM, the interaction is mediated by the dynamical charge susceptibility of the Fermi gas, Il{q — l,ujp — oji) 
(see Fig. ^(b)), which is generally a complex function of momentum and frequency. This gives rise to two effects 
which both are relevant to our analysis: 

(i) n(g — ljUp — uji) has a branch cut, which gives rise to a nonzero renormalization of the spin wave dispersion at 
T = 0, and 

(ii) the static Il{q — 1,0), relevant to thermal corrections to the dispersion at T < JiS, scales differently from Ji 
and this changes the scale of the self-energy in the DEFM relative to a classically equivalent HEM. In particular, at 
smallpi?, Ji oc pp, while Il{q — l) oc pp at q,ltyp <C pf and Il{q — l) cx p'^p/\q — l\'^ for \q — l\ ^ pf- Thermal corrections 
to spin wave dispersion in the DEFM are therefore enhanced relative to the HEM. Furthermore, at small doping x, 
when typical q and / exceed pf, the momentum dependence of n(g — ^) causes thermal self energy correction in the 
DEFM to have a different functional form from that in the HEM. 

We now proceed with the calculations. Assembling all three 0{1/S'^) contributions to the self-energy and splitting 
it into quantum and thermal pieces, we obtain after some algebra that all unwanted terms are canceled out and the 
Goldstone mode at q = survives, as required. The resulting self-energy is given by 



I](2)(g,c.) « S(2)(g) = 42^^(9) + (g) 



(47) 



where 




zt 1 



pi 



■p 



Ip+li 
Ip - 11 



Ip 7fc+; 
Ip - 11 



(48) 



and 
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FIG. 4. Dispersion of spin waves in a cubic DEFM with electron doping x = 0.7 and spin S = 3/2, in units of the electron 
bandwidth 2zt = 1.0. Solid line — dispersion of classically equivalent Heisenberg model. Dashed lines — dispersion of DEFM 
including leading quantum and thermal corrections for (top to bottom) T = 0, T = 2D/3 and T = D where D = 12715 is the 
bandwidth of spin wave excitations. 



zt 1 



{ip - lp+q){lp+q-l - Ip+q) , [ip " lp+l)[lP + I 



Ip+l)' 



Ip - Ip+q-l 



Ip - lp+l~q 



(49) 



Here np{k) = npiik) and nsiq) = nBi^iq)) are Fermi and Bose distribution functions, respectively. The frequency 
dependence of the self-energy is relevant for the computations of the bosonic damping which is of order 1/5*^ (see 
below), but can be neglected in the computations to order l/S*^ as typical bosonic uj are of order uJsf{q) oc t/S and are 
small by 1/5' compared to typical fermionic energies which are of order t. In this situation, the full bosonic dispersion 
is simply 



UJsiviq) = zJiS{l - 7g) + S^(q). 



(50) 



We will also neglect the temperature dependence of the Fermi functions as we are interested in temperatures of order 
the spin wave bandwidth T ^ t/ S <^t. This makes possible the simple separation of the self energy into thermal and 
quantum pieces as written in Eq. (47). We now analyze quantum and thermal corrections to spin wave dispersion. 
We begin with the case of zero temperature. 



A. Quantum corrections at T = 

1. spin wave dispersion 

We see from (^) that at T = 0, the interaction with the charge degrees of freedom leads to an overall reduction in 
the spin wave bandwidth (the first term in (48)), together with a modification in the form of the dispersion provided 
by the second term in in ([48|). This second term 



pi 



(51) 



is either positive or negative throughout the Brillouin zone, depending on doping, and has the symmetry AS^^^ [q) - 
AE(^^(7r — q) which bare spin wave dispersion does not possess. Near the center of the Brillouin zone it behaves as 



AS(^H<Z«0)~ 252^(a;)<z2 



(52) 
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FIG. 5. Numerical results for the imaginary part of D{q, uj) of a 2D double exchange ferromagnet with 5* = 3/2 and doping 
of 0.18 holes per site, with energy in units of the hopping integral t. The plots show the increasing damping of spin waves along 
the zone diagonal. Sharp structure in the data is an artifact of finite numerical resolution and has been partially removed by 
convolution with a Gaussian of width O.Qlt. 



where 

^(^) = ^E^^[l-"^(^2)]nF(A;i) (53) 
N'^ ^ 71 - 72 

is a doping dependent constant. By symmetry, the behaviour of AI]'-^-'((7) near the zone corner q k, i: must have 
exactly the same form. Along the zone diagonal 

AI](2)(g) = ^/(x)[l-cos(2g)] (54) 

This form is vepjij similar to the form of correction to Heisenberg dispersion found in numerical studies of the DE 
model on a ringla. 

For a small density of electrons {x < 0.31), I{x) > in both 2D and 3D, and corrections to dispersion lead to 
a "softening" of modes at the zone boundary relative to those the zone center. The same is true for small density 
of holes 1- X < 0.07 in 21? and 1 - a; < 0.06 in 3D. For intermediate densities 0.31(7) < a; < 0.92(7) in 2D, and 
0.31(7) < X < 0.94(2) in 3D, I{x) is negative, i. e. quantum effects cause a relatiice softening of spin wave modes at 
q = 0. These results are in good agreement with earlier numerical studies studiesEHI. We discuss the consequences of 



the non-monotonic doping dependence of I{x) in more detail later in Sec. IV when we consider order from disorder 
effects. 

We note in passing that if we formally extend our 1/S results to arbitrary S, the spin wave dispersion becomes 
unstable for S < Scr where in 3D Scr — (1/2A^) J2ki^ + 7fc)/(l — 7fc) ~ 1- This opens up a possibility tiat for small 
S (e.g., S = 1/2), the ground state may not be a ferromagnet, as suggested by some numerical studiesli3. However, 
the extension to small S is beyond the scope of the present paper. 



2. Damping of spin waves 

While spin waves are exact eigenstates of the Heisenberg Hamiltonian (||), in the DEFM only the q = (Goldstone) 
mode is an eigenstate. All other spin wave modes in the DEFM have a finite lifetime, even at zero temperature. 
Physically, the possibility for a spin wave to decay at T = is again related to the presence of Fermions. A spin 
wave with energy ajswil) can give energy to a particle-hole pair, and decay into a lower energy spin wave — the 
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process described by the diagram Fig. 113. As a consequence T,^ '(q,fl) contains an imaginary part which describes 
the damping of a spin wave. 

To calculate damping effects it is necessary to keep the bosonic energies i7, lo^ together with the fermionic in 

the denominator of S^^)(g, 51). This means that the calculations are formally less controlled, as we include terms one 
order smaller in 1/5*. To work rigorously at this order it would be necessary to extend the (inverse) transformation 
Eg. (pO|-p^) to 0{1/S^) — a very involved task. However the physical origin of the damping is clear, and we verified 
that the imaginary part of the self-energy to f{l/ S^) is fully captured by just keeping il in denominators of two-loop 
self-energy diagrams. We then obtain 

^T^oiq, ^) = Q(f^) (4^)2 ]^ E + - 2efe+,)' npik + q')[l - nFik)]5 {{n - u;,^,,) - {e^+g' - e^)) (55) 

where we have grouped together terms of 0(0) and terms of 0(1/5") within the argument of the delta function. 

Without S"((7, fl), the imaginary part of the bosonic propagator D{q, 57) is just a delta function dispersing with the 
spin wave energy ^^111(9), i- e. 

D"{q,n) = -nSin-uJsy,iq)) (56) 

However because of the non-zero Y}"{q,fl), the simple delta function peak ( |56| ) is replaced by a broader dispersing 
feature with is approximately Lorentzian in shape with the maximum at uJsw{q)- 

We can obtain an estimate of the width of D"{q, 57) by calculating the damping of spin waves on the mass shell 

r? = E (^'^ ~ ^F{k + q')[i - nF{k)]S {(ujsw{q) ~ i^sw{q - q')) - {fk+q' - tk)) ■ (57) 

kq' 

where we have used the delta function on energy to eliminate terms of subleading order in l/S* from the vertex. The 
damping of spin waves must vanish for 5 ^ 0, in order for the Goldstone theorum to be satisfied. For q <^ pp we 
find that F^ vanishes as u!sw{q) q"^^^ ^ q"^^^. A factor ujsw{q) comes from the 2D integration over d\k — kp\ and d\q'\, 
another factor q^ comes from the vertex, and the remaining factor q"^^^ comes from the fact that \q'\ < q. nad hence 
d'^g' ~ q''-^^d\q'\. In 2D, this yields F, oc q^ and in 3D, F^ cx q^. This fully agrees with earher estimatescllE^I. This 
strong dependence on q survives wp to q ^ pp- At larger q, and pp <^ 1, the damping of spin waves saturates at a 
constant value F^ oc (t/S^). At pp '^-^ 1, F^ at 9 > pp is in general a rather complicated function of momentum. 

In Fig. 1^ we plot the imaginary part of D(q^ Q) obtained by numerically evaluating (|5^) for a set of momenta on 
the zone diagonal for a DEFM on a 2D square lattice, for S = 3/2 and electron density x = 0.82. From the plot, 
we see that zone corner modes are clearly broadened relative to those in the zone center. Evaluating F^ througho«jt 
the BZ for a range of dopings in 2D we obtain a functional form of damping similar to that recently reported inll3. 
However our estimate on the upper bound for values of damping is about 25% smaller than that given ir£^. 

The results for the 3D cubic case are similar, and best summarized by plotting the F^ along all major symmetry 
directions. In Fig. ^ we plot the damping for a cubic system with electron doping x = 0.7 and spin S = 3/2, in units of 
the electron bandwidth 2zt = 1.0. The damping of spin waves is small at low q, but strongly momentum dependent. 
Damping is large for q approaching the zone edge, and has a maximum value of about 10% of the spin wave dispersion 
at the zone corner. Stationary points of the damping (maxima, minima and points of inflection) occur at symmetry 
points of the the BZ or where q crosses the Fermi surface. 

The large absolute value of the damping at large q in the DEFM and its strong momentum dependence is consistent 
with the experimental behavior seen in Neutron scattering experiments on the Manganites, wheiie damping is large, 
highly momentum dependent, and can rise to ^ 10% of spin wave dispersion at the zone cornerO. 



B. Thermal corrections at T > 

We now discuss the form of the correction to the spin wave dispersion at finite T. Continuing the spirit of comparison 
with the low T expansion for a Heisenberg FM, we assume that T < JiS <^ t. In a nearest-neighbor Heisenberg FM 
the thermal renormalization of the spin wave dispersion to first order in l/S is 

where0 
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FIG. 6. Damping of the spin wave on the mass shell as a function of momentum throughout the BZ for a cubic DEFM with 
electron doping x = 0.7 and spin 5* = 3/2, in units of the electron bandwidth 2zt = 1.0. The volume enclosed by the Fermi 
surface is denoted by shading. 



B 

At the lowest temperatures, 



We see that the dispersion preserves its T = form, and that thermal fluctuations only reduce the overall scale of 
dispersion by an amount which at low T is proportional to T^^^. For the DEFM we find below that at small x, the 

(2) (2) 

form of Ej, (g) is totally different from that in HFM. For finite x, we find that {q) can be crudely approximated 
by the same functional formulas (^0|), but the overall factor has completely different doping dependence from the 
effective exchange coupling Ji. 

We now proceed with the analysis of Eq. (^). In general, at finite T there are three typical momenta in the 
problem, the spin wave momentum q, the Fermi momentum pp, and a thermal length scale Ityp ~ (T/ JiSy^^ (in 
units where the lattice constant a = 1). Let us consider first the case when T is very low and Ityp ^ (T/JiS")^/^ are 
smaller than pp. This case mimics the behavior at intermediate electron densities. Expanding in (^) in I we obtain 
after a simple algebra: 

^ril) = ^Y.''B{l)l'Y.^F{p)lp^{p.q) (61) 



4^2 



where 



= 1 + 7 - — V ^"^^(P' + + sin^(p^ - qi)+sin^{pi) 

where i takes the values x, y and z. If ^'(p, q) was just proportional to 1 — 7^, the thermal self-energy would have the 
same functional form as in HFM. We see, however, that ^(p, q) is more complex and depends on both variables. 
The form of ^'(p, q) is particularly simple in ID. Here we obtain from Eq. (^2|) 



^{p, g) = (1 + cosq) '^"'(^ + '^^ + '^"'^^ - '^^ + 



cosp(cosp — cos(p + q)) cospicosp — cos{p — q)) 
l + 2sin^p + cos(7 

(l-cosg)- . (63) 

1 — cos g — 2 sm p 
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We see that at small q and finite p, 'i>{p,q) scales as as it indeed should. Substituting the result for '^{p,q) into 
( |6l| ) we find after integrating over p 

+ I \ 

- cos q)N{q) (64) 




where 



At g ^ 0, 



sinpF / , , 1 , sinpF + sin(7/2 , 

TT \ Zsmppsmq/Z \smpF — smq/ Z\ 



At the zone boundary, g = tt, 



N{q) ^ - f sinpF ) . (66) 

TT \smpF 



sinpF / , 1 , l + sinpF\ .„„^ 
7V(7r) = -1 + ^ fog-^ , . (67) 

TT \ ZsmpF 1 — smpF, 



At smallpF, A''(7r) = (I/Stt) sin + We see that N{q) > 0, and hence (q) is negative, i. e. thermal fluctuations 
reduce the spin wave energy. This agrees with the spin wave result for the HFM. 

In Fig. we plot the form of thermal corrections to spin wave dispersion for a DEFM in ID with pF ~ 0.4, together 
with a correction proportional to 1 — 7^ having the same prefactor at small q. We clearly see two effects — the 
suppression of temperature corrections at large q relative to small q, and a logarithmic singularity coming from the 
perfect nesting of the Fermi surface in ID. 

The q^ behavior at the smallest q exists in all dimensions as we have explicitly verified. The explicit forms are, 
however, rather complex and we refrain from presenting them. It is essential that for all I?, the prefactor scales as 
1/p^ for g <C p <C 1. In 3D, the q^/p^ form of 5'(p, q) yields for g, ^ <C < 1, 

4^^(.) - -^^^E-BUr « q'p.T^^' (68) 

In 2D, we obtained 
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q = (0,0, It/2) o 
q = ( , , JI) + 
q = (IC/I, Il/2,Jl/2) □ 




FIG. 8. Ratio of thermal corrections to spin wave dispersion in a DEFM to the thermal corrections of a classically 
equivalent Heisenberg FM a{q,x,T) = Y,^'^''^' {q) /Y.^^'^ (q) for momenta g = (0,0,7r/2), q = (0, 0, tt) and q = (7r/2, 7r/2, 7r/2) 
as a function of doping x and reduced temperature T/JiS. Plot shows that data for a wide range of dopings, momenta and 
reduced temperatures approximately collapse onto a "universal" surface. 



^'T\l) = -^,^T.-Biir^,'T^ (69) 

Since the overall scale of spin wave dispersion in the DEFM at low doping is proportional to x, (which scales as 
in 2D and in 3D), we see that the overall scale of thermal corrections in the DEFM is enhanced by a factor l/pj^ 
relative to a Heisenberg model with the effective FM coupling Ji. 

Numerical results for the spin wave dispersion in a cubic system with an electron doping of a; = 0.7, including thermal 
corrections, are shown in Fig. |^. Thermal corrections are generally small compared with the quantum corrections at 
T = 0, the two effects becoming comparable only for T ~ 6Ji5. 

It is helpful to divide thermal self energy corrections in the DEFM by the thermal corrections for a Heisenberg 
model with the same spin stiffness JiS, i. e. to consider 

We plotted this ratio against reduced temperature T/JiS for several q and and a range of dopings 0.1 < a; < 0.9 in 
Fig. ^. We see that for a wide range of dopings and reduced temperatures the data for different q collapse onto a 
single "universal" surface. This means that a remains roughly constant as a function of T/ JiS, and also depends only 
weakly on on q, which implies that for the experimentally relevant intermediate dopings, the thermal self-energy in a 
DEFM roughly mimics that in a HFM. Still, however, the overall amplitudes of the corrections are very different for 
the reasons given above. We found that the ratio of the two is approximately constant at a sa 5 for 0.25 < x < 0.75, 
falling away to about a = 2-3 for x = 0.1 {x = 0.9). These numbers should be compared to enhancement of corrections 
by-a factor 2-3 relative to a nearest neighbour HFM observed in Lao.ssPbo.isMnOsES, for which x = 0.85. The authors 
oiE3 interpreted the enhancement of temperature corrections in Lao.ssPbo.isMnOs in terms of effective non-nearest 
neighbour couplings — we note that these are dynamically generated by the form of interaction between spin waves 
in the DEFM — see Fig. ||(b). At the smallest dopings a acquires a strong dependence on both q and T/JiS, as 
discussed below. We also found that for Z? = 2 and D — 1, the deviations from a « const are more prominent for 
the same intermediate x. In particular, in ID a simple experimentation with trigonometry yields N{tt) < N{0). This 
implies that compared to the HFM, the strength of thermal fluctuations in DEFM is reduced near q — tt. 

We now consider analytically what happens if we lift the restrictions that both q and Ityp T / JiS are smaller than 
Pf- As we still focus on low T < JiS, pp should also be small which in turn implies that electron density a; <C 1. We 
verified that the results similar to the ones below also hold for small hole density 1 — a: <C 1 . 
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FIG. 9. Thermal corrections to spin wave dispersion in a cubic DEFM with x = 0.67 for q = (7r/2, 7r/2, 7r/2) (points) 
compared with those for a classically equivalent Heisenberg FM (solid line). While the functional form of the corrections is 
similar, their magnitudes are very different. The inset shows the crossover from y^/^ behavior at low T to T^/^ behavior at 
T > JiS. 



Let us first consider what happens when q exceeds pp- In this limit, we can approximate ^'(p, g) by '^(O^q). 
Substituting the resuh into the expression for the self-energy, we obtain, in any D 



(71) 



where 



In ID, (171^ reduces to 



4'^ (9) 



2 1 



l2q 



-(1 + 7,) 



tx\ 



Vb(0/^(1 + cosg) 



(72) 



(73) 



This result coincides with Eq. ( p4[ ) as ior pp <C 1 and q ^ pp, N{q) from ( |65[ ) reduces to N{q) — x{l + cos q) / {1 ~ cosq) 
where x = pp/n. 



In 3D, Eqn ^) yields S^^^(g) cx p^T^/^ 



In 2D, we have H^^q) oc PpT^. Comparing these results with Eq. ( |6§| ) 

and (|69D, we see that the at q ^ pp, the self-energy is reduced by p'^/q^, relative to corrections at the zone center. 

This is precisely the result we anticipated on general grounds as when q exceeds pp, the effective boson-boson vertex 

is reduced due to the reduction of the charge susceptibility which mediates boson-boson interaction. Note that this 

reduction eliminates the quadratic dependence on q, i. e. the thermal self-energy becomes flat at q^ pp. 

Curiously enough, in 3D, at these low dopings, P{q) is negative for all q along zone diagonal {P{q) = — (1/3)(1 — 

cosg)), i. e. the thermal self-energy changes sign compared to Eq. (^). In 2D, P{q) vanishes along zone diagonal, 

and in ID, P{q) = 1-1- cosg is positive shown above (see Eq. ([73|). 

The reduction of the bosonic self-energy due to the reduction of the charge susceptibility at large momenta can 

(2) . 

also be detected by analyzing the form of when q is smaller than pp, but Ityp exceeds pp. This case is even more 

(2) . . 

instructive than the large q and small Ityp limit, as can be obtained analytically for arbitrary ratio of Ityp and 
pp. Indeed, expanding in all three momenta in (p9[) we obtained in 3D 



tq^pp 



1 



1447r2S'2 N ^ 



(- 

\2pp 



(74) 
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FIG. 10. Thermal corrections to spin wave dispersion in a DEFM with x = 0.67 normahzed to those in an equivalent HFM 
for q = (0, 0, 7r/2), q = (0, 0, 7r),g = {tv/2, tv/2, 7v/2) and q = (vr, tt, tt). 



The scaling function (y) is given by 



(75) 



For y < 1, Zmin — 0, for y > 1, z„iin — \/ (y^ — ^)/y'^- Evaluating the integral analytically in the limiting cases, we 
find $3£)(y ^ 0) = 1, $3D(y = 1) = 0, <f>3D(y > 1) ~ -1/j/^- The y = limit reproduces Eq. (||). For large y, 
i. e. for g < < {T/JiSY'^, we have S^^^(g) = Cq^p%{T / JiSf^ where C > 0. The scaling function $3d(?;) is 
shown in Fig. (|ll|). The sign change of <i>3D(y) between small and large y collaborates our earlier observation that 
in 3D, the self-energy at the smallest T and at T/JiS > pp have different signs. When q exceeds pp and become 
comparable to (T/JiS*)^/^, Eq. ( [zil) becomes invalid, and Y]^^\q) smoothly interpolates into T,^^\q) (xp'pT^/'^, as in 
Eq. 0. 

A very similar behavior holds in 2D case. Here we obtained 



4'^ (9) 



tq^ 



487r52 N 



\2pF 



(76) 



where 



*2d(2;) = 1, y<i 

Try^ Jo 



Vl- (2/sin0)2(l + 2y2 + 4(?;cos0)2), y>l 



(77) 



l/y. In the limit ?/ 3> 1, we found analytically <i>2D(j/) = ci/y'^ with a > 0. 

Nl/2 



where sin (l},nax 

Again, when q exceeds and becomes comparable to : {T/ J*S)^^'^, Eq. ( |7^ ) becomes invalid, and E^-'(g) smoothly 

.(2) 



interpolates into E^^^(g) cx pf^T^/^ Eq.(|7l]) 



Finally, in ID, 



£<?>(,) 



2p7 



(78) 



where 
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FIG. 11. Scaling function governing thermal corrections to spin wave dispersion in the DEFM in ID — Eq. (|79|), 2D — 
Eq. ^ and 3D — Eq. (Q. 



In the two limits, $i£)(0) — 1, and ^ 1) = before, when g exceeds pp and becomes comparable to 

{T/J*SY/^, Eq. @ smoothly interpolates into Eq. @. 

We see that in any D, the Heisenberg form of the correction to the spin wave dispersion is reproduced only at 
small {T/JiSy^^ ^ pp, when the scaling function can be approximated by a constant. At larger temperatures 
ytyp (T/JiSpj^y^^ become large, and the variation of cannot be neglected. The only difference between 

different D is that <&3D(y) changes sign at y = 0{1) and becomes negative at large y, while <i>2D and remain 
positive for all y. 

Unfortunately, the new physics we just described is confined to very small x'm D=3. This is simply because in 3D 
Pf = {Gn'^xy^^ , and to make pp small, one needs a very small x. For such small x, our numerical analysis is not very 

accurate, but we were able to confirm that (g) changes sign as T increases. We also clearly detected the downturn 

renormalization of |E^'^(g)| with increasing T at x (or 1 — x) ~ 0.1. 



IV. ORDER FROM DISORDER IN THE DOUBLE EXCHANGE MODEL 



We now apply the results of the previous section to study the order from disorder phenomenon in the double exchange 
model. We recall that if double exchange ferromagnetism Ji is put into competition with direct antiferromagnetic 
superexchange J2 between core spins, then classically, ferromagnetic state becomes unstable at Ji < J2, but the 
intermediate configuration is not specified — the only requirement on classical spins is that neighboring spins are 
misaligned by an angle 6: 

cose = Ji/J2- (80) 

This condition allows an infinite set of intermediate states. The two limits of this set are canted and spiral spin 
configurations, as discussed in the introduction. 

We first demonstrate that at a classical level, the spin wave excitation spectrum of either canted or spiral state indeed 
contains a zero-energy branch of excitations. Consider for definiteness a canted phase. The bosonic Hamiltonian for 
a two sublattice canted structure described by the full TYq + Tii + TL2 at vanishing t/ Jh and J2/ Jh can be obtained 
directly within our expansion scheme. Working to 0{1/ S) and introducing two quantization axes misaligned by 0/2, 
we obtain 

He = nf^ + nf^ +nf^ + ... (si) 
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where 



7^« = 



E 



J2S^ COS 61 



(82) 



sm — 
25 2 



/L/b, - al) + h. 



-J2S\ — sin 



(2) 



E 

+J2S 



45^°^ 2 



} 



2aj6, ) + /i. c. 



(^dlbj + 6]a,) + sin^ - (a,6j + fo^gt 



(83) 



(84) 



where /a denotes a local "up spin" electron on sublattice A, b a, Holstein-Primakoff boson on sublattice B, etc., 
and the sum over (ij) is restricted to pairs of neighbouring sites with site i on the A sublattice and site j on the B 
sublattice. We note that the effective Hamiltonian for the antiferromagnetic Neel state is obtained by setting — it 
above. 

To zeroth order in 1/5* the ground state energy of this configuration is 



Eq = 2Nz 



J2S cos 9 — t cos ■ 



(85) 



Minimizing this with respect to 9 gives the de Gennes result (Eq. |o|) for the canting angle. 

The spin wave dispersion to leading order in 1/S is obtained in the same way as for a ferromagnet, by evaluating 
the leading corrections to an (initially dispersionless) bosonic propagator. This spin wave "spectrum" is doubly 
degenerate, and the lattice momentum q is restricted to the reduced Brillouin appropriate to the two sublattice spin 
configuration. The evaluation of the self-energy to order l/S* is technically a bit more involved than for a ferromagnet 

as in addition to the one-loop self-energy from Tif^ we also need to compute the second-order contribution from Tig^''. 
The last term vanishes at zero momentum after one substitutes the canting angle, as it indeed should, but it does 
contribute at second order. We found that the two contributions exactly cancel each other, i. e. to order 1/S, 



(86) 



The vanishing of the leading order spin wave dispersion is clearly linked to the degeneracy of the classical groundstate, 
and, as we already mentioned in the introduction, has a simple physical interpretation — at the classical level we 
may take any spin on the B sublattice and rotate it about the direction of magnetization of the A sublattice without 
altering the angle between it and the neighbouring spins. As doesjipt change, this perturbation does not cost any 
energy, and by extension the spin wave dispersion at 0{1/S) is zeroB'El. 

We now proceed to the computations beyond the leading order in 1/S. Our key goal is to understand which 
configuration of the classically degenerate set is selected by either quantum or thermal fluctuations. There are two 
ways to attack this problem. First, we may continue with the analysis of the canted phase and perform computations 
to next order in 1/5. This is possible but requires a lot of efforts as at each order there is a competition between the 
contributions from from the term quartic and bosons and the term linear in bosons. A second route which we choose 
is to analyze the momentum of the instability of the ferromagnetic configuration. To 0(1/5), the excitation spectrum 
of a ferromagnet in the presence of the antiferromagnetic exchange is given by Eq. (R) which vanishes identically 
for Ji — J2, in agreement with the infinite degeneracy of the classical groundstate for Ji < J2. We already know, 
however, that the excitation spectrum of the DE model is different from that of a ferromagnet and, depending on the 
density of electrons, is softened either at the zone boundary or at g = 0. 

A way to understand what this implies is to analyze what form of the spin wave dispersion one generally expects 
at a point where a ferromagnet becomes unstable towards either canted or spiral phases, if there were no degeneracy. 
The simplest way to eliminate a degeneracy is to formally add an exchange interaction between second neighbors 



H3 — J3 SiSj 
This extra interaction modifies the dispersion Eq. (^ to give 



(87) 
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LUswiq) = zS{Ji - J2)(l - 7q) - zSJj,{l - J2q) 



For positive J3 (antiferromagnetic next-nearest neighbour exchange), the J3 term favors maximum misahgnment 
of the next-nearest neighbour spins. This condition is best satisfied by the spiral state. For negative J3, the extra 
couphng favors parallel orientation of next-nearest neighbour spins. This is a property of the canted phase. Analyzing 
the dispersion in the presence of J3 we immediately see that for positive J3, the spectrum first becomes unstable at 
g = 0, at Ji = J2 + 4J3, while for negative J3, the excitation spectrum first softens at Ji — J2 aX q — (7r,7r,7r). 
Performing the classical spin wave analysis for the full DE Hamiltonian with extra J2 and J3 terms slightly below the 
instability we find that for positive J3, the spin wave spectrum has three zero modes: at g = and an incommensurate 
momenta q = ±(go7 ^Oi 90) where = 6. This is exactly what we expect in a spiral state. For negative J3, the soft 
modes remain at g = and at (tt, tt, tt). This is what we expect in a two-sublattice canted phase. 

The above analysis shows that the type of the ground state configuration below the instability can be obtained by 
just analyzing the spin wave spectrum at the instability point of the ferromagnetic phase. Using the results from the 
previous section we find that at T = 0, the instability occurs when 



(89) 



Cos^o{q) = zS(Ji - J2)[l - 7(9)] + AS('n'?,0) = 

where, as shown, Ji = t/2S'^ with t = tN~^ Ukjk, and the relevant self-energy correction is presented in Eq. (|4J 
As discussed in Section this can be split into two of two pieces. One scales as 1 — 7(g) and just renormalize Ji. 
The second piece AE^^^ [q) accounts for the renormalization of the form of the dispersion. It reads 



zt 



-y 



7f - 7i- 



-ni(l - ^2) 



(90) 



The first instability of the DEFM will occur for some value of J2 sufficient to cancel the bare spin dispersion and 
this quantum correction. The value of critical J2 depends on the form of AS*-^-* as a function of g and doping x. We 
therefore consider the ratio 



A? 



AI](^)(g) 
1-7(9) 



(91) 



as a function of doping. The wavevector g = g* for which has its minimum value defines the momentum of the 
instability. 

Near g = 0, we obtain using ( p2|) and (|5 



t I{x)q^ 



252 g2/z 25' 



while at the zone corner 



A^. 



t I{x)q^ 
252 2-(7r-g)2/2 







(92) 



(93) 



discussed in Sec |I] 
and 1 



Whether the first instability occurs in the zone center or the zone corner therefore depends on the sign of I{x). As we 

^ ^ I{x) is positive for x < 0.31 in both 2D and 3D, and for small density of holes l — x< 0.07 in 2D 

X < 0.06 in 3-D. For these dopings, the intermediate configuration is then a canted phase. For intermediate 
densities 0.31(7) < a; < 0.92(7) in 2D, and 0.31(7) < x < 0.94(2) in 3-D, I{x) is negative, i. e. quantum effects cause 
a relative softening of spin wave modes at g = 0, and the intermediate configuration is a spiral phase. These results 
are in good agreement with earlier numerical studies studiescj. 

At finite temperature, the range of dopings where the canted phase is selected increases with T as for realistic x, 
classical fluctuations primarily soften the dispersion near g = and hence favor the spiral state, at least at small T. 
At higher T, and in -D = 3, thermal fluctuations may in principle, cause a re-entrant transition into a spiral phase 
as when T becomes larger than JnSp^, the sign of T,^^\q) changes, and thermal fluctuations now favor the canted 
phase. In practice, however, we found that this sign change occurs at unrealistically large T for which our low T 
analysis is inapplicable. 

We conclude this section with few words of caution. First, in our analysis we assumed that there are no nontrivial 
transitions inside the intermediate phase, i. e. the configuration which becomes the ground state immediately after 
ferromagnetic instability yields to another configuration at some distance away from the instability. Such a possibility 
exists in general, Jiut to the best of our knowledge, there is no known example of such behavior. Second, we neglected 
phase separation^. To study this possibility in our approach, one has to analyze the sign of the longitudinal suscep- 
tibility in, e.g. the canted phase. If it is negative, then the system is unstable towards phase separationcJ. These 
calculations are currently under way. 
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V. CONCLUSIONS 



We have performed a novel large S expansion for the Kondo model on a magnetically ordered lattice, and presented a 
simple, physically transparent and controlled calculation scheme for evaluating corrections to magnetic and electronic 
properties of the model to 0{1/S'^). Calculations are particularly straightforward in the limit where the onsite 
(ferromagnetic) Kondo coupling Jh greatly exceeds the electron bandwidth zt (the double exchange model). In this 
limit we have used our expansion scheme to study quantum and thermal corrections to the spin wave spectrum of 
the double exchange ferromagnet. We argue that a double exchange ferromagnet is equivalent to a Heisenberg model 
with effective coupling Ji oc t/S'^ only at a quasi-classical level — both quantum and thermal corrections to the 
spin wave spectrum of a double exchange ferromagnet differ from any effective Heisenberg model because its spin 
excitations interact only indirectly, through the exchange of charge fluctuations. We demonstrated that the effective 
interaction between spin waves is mediated by Pauli susceptibility of electrons which has a much larger amplitude 
than Ji and also has its own non-analytic momentum and frequency dependence. We argued that at T = 0, the 
nonanalyticity in the frequency dependence of the effective vertex gives rise to a finite self-energy correction which 
changes the momentum dependence of the spin wave dispersion compared to that in a Heisenberg ferromagnet. The 
form of the full dispersion which we found is similar to that found in numerical studies of the DE modelta. For a small 
density of electrons, or small density of holes, the corrections to dispersion lead to a "softening" of modes at the zone 
boundary relative to those the zone center. For intermediate densities, quantum effects cause .a, relative softening of 
spin wave modes at g = 0. These results are in good agreement with earlier numerical studies c3. 

The softeaing of zone boundary spin waves has been observed in neutron scattering experiments on the CMR 



Manganitesca at x ^ Aif-i 1^ been suggested that this_softening may be due to deviations £mm Heisenberg 



behavior in the DEF the influence of optical phononsE3, or due to orbital degrees of freedomcZl. Our analysis 

argues against the first possibility as for x ~ 0.7, our zero temperature theory predicts a relative hardening rather 
than softening of the dispersion at (7r,7r,7r). We note however that the leading corrections to spin wave dispersion at 



finite t/Jn do lead to a softening of zone boundary mode£3. 

We also found that the process by which spin waves decay into lower energy spin excitations dressed with particle- 
hole pairs leads to a finite spin wave lifetime at T = 0. The spin wave lifetime is very long at low q, scaling as in 
3D for > g 0, but increases rapidly with increasing q. The large absolute value of the damping at large qdn the 
DEFM and its strong momentum dependence is consistent with the experimental behavior in the ManganitesEa. 

We then analyzed in detail the form of the temperature corrections to the spin wavespin wave dispersion. We argued 
that for experimentally relevant densities, the corrections roughly have the same functional form as in a Heisenberg 
ferromagnet, but the overall scale is much larger. At low density (small pf), for which an analytic treatment is 
possible, thermal corrections in DEFM are enhanced relative to those in a HFM with the same effective exchange 
coupling Ji by a factor proportional to l/p%- At very small pi?, we also found that not only the amplitude but also the 
functional form of the thermal self-energy in DEFM is qualitatively different from that in a HFM. In 3D, the thermal 
self-energy in DEFM at low doping even changes sign relative to that in a HFM. The enhancement oLtemperature the 
overall scale corrections relative to a HFM has been observed experimentally in Lao.SSSro.lSMnOstj, and we believe 
that our finite temperature results will be useful for the interpretation of the neutron data on Manganites. 

We also considered the effect of a direct superexchange antiferromagnetic interaction between core spins. We find 
that the competition between ferromagnetic double exchange and an antiferromagnetic superexchange provides a 
new example of an "order from disorder" phenomenon — when the two interactions are of comparable strength, 
an intermediate spin configuration (either a canted or a spiral state) is selected only by quantum and/or thermal 
fiuctuations. We discussed which configurations are selected at various dopings. 

The issue left for further study,is the possibility of a phase separation in the intermediate regime into ferromagnetic 
and antiferromagnetic regionsQca, and also possible stripe formation at low hole doping. The analysis of these issues 
is clearly called for. 

Acknowledgments. 

It is our pleasure to acknowledge helpful conversation with S. W. Cheong, D. Golosov, R. Joynt, M. Kagan, T. A. 
Kaplan, G. KhaliuUin, S. D. Mahanti, A. Millis, E. Miiller-Hartmann, M. Norman, O. Tchernyshyov, N. B. Perkins, 
N. M. Plakida and M. Rzchowski. This work was supported under NSF grant DMR-9632527 and the visitors program 
of MPI-PKS (N. S.), and by NSF grant DMR-9979749 (A. V. C). 



20 



^ C. Zener, Phys. Rev. 82, 403 (1951). 

^ P. W. Anderson and H. Hasegawa Phys. Rev. 100, 675 (1955). 

3 P. Majumdar and P. Littlewood, Phys. Rev. Lett. 81, 1314 (1998). 

* J. Riera, K. Hallberg and E. Dagotto, Lett. 79, 713 (1997). 
^ P.-G. de Gennes, Phys. Rev. 118, 141 (1960). 

D. Golosov, Phys. Rev. B 58, 8617 (1998). 

''A. B. Harris, C. KaUin and A. J. Berhnsky, Phys. Rev. B 45, 2899 (1992); S. Sachdev, Phys. Rev. B 45, 12377 (1992); A. 
V. Chubukov, Phys. Rev. Lett. 69, 832 (1992). 

* Nic Shannon and Andrey Chubukov, J. Phys. Condens. Matt, 14 L235 (2002). 
" Nic Shannon, J. Phys. Condens. Matt, 13 6731 (2001). 

This is indeed true for the computations of the transverse spin correlator between different sites. For a Kondo problem at a 
given site, higher-order terms in the Holstein-Primakoff expansion for a composite spin are indeed relevant as they preserve 
the correct spin algebra, and are incorporated into our transformation in Sec. |l^. 
" K. Kubo and N. Ohata, J. Phys. Soc. Jpn. 33, 21 (1972). 

E. L. Nagaev Phys. Rev. B 58, 827 (1998); 

"Physics of Magnetic Semiconductors", Moscow, Mir Publ (1979). 
" N. Furukawa, J. Phys. Soc. Jpn. 63, 3214 (1994). 

A. J. Millis, P. B. Littlewood and B. L Schraimann Phys. Rev. Lett. 74, 5144 (1995). 
N. B. Perkins and N.M. Plakida, Theor. Math. Phys., 120 (1999) 1182. 

D. Golosov, Phys. Rev. Lett. 84, 3974 (2000). 

see e.g., Akhiezer et al "Spin waves", Wiley, NY (1968). 

T. A. Kaplan and S. D. Mahanti, J. Phys. Condens. Matt. 9 L291 (1997). 

J. Zang et al, J. Phys. Condens. Matt. 9, L157 (1997) and references therein. 

P. Wurth and E. Miiller-Hartmann, Europhys. J. B 5, 403 (1998); D. Golosov, J. Appl. Phys. 87, 5804 (2000). 
2^ V. Yu Irkhin and M. 1. Katsnel'son JETP 88, 522 (1985). 

22 L. Vasiliu-Doloc et. al. Phys. Rev. B 58, 14913 (1998). 

23 A. Moreo, S. Yunoki and E. Dagotto, Science 283, 2034 (1999). 

B. L Shraiman and E. Siggia, Phys. Rev. B 46, 8305 (1992); A. Chubukov and K. Musaelian, Phys. Rev. B 51 12605 (1995). 
2® H. Y. Hwang et. al. , Phys. Rev. Lett. 80, 1316 (1998), Tapan Chatterji, private communication. 

2^ N. Furukawa, J. Phys. Soc. Jpn. 68, 2522 (1999). 

2'' G. Khahullin and R. Kilian, Phys. Rev. B 61, 3494 (2000). 

2* F. Mancini, N. B. Perkins and N. M. Plakida Phys. Lett. A 284, 286 (2001); Nic Shannon, in preparation. 
23 M. Yu. Kagan, D. L Khomskii and M. V. Mostovoy Eur. Phys. J. B 12 212 (1999). 



21 



